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We propose a new locally smeared operator product expansion to decompose nonlocal operators 
in terms of a basis of smeared operators. The smeared operator product expansion formally con¬ 
nects nonperturbative matrix elements determined numerically using lattice field theory to matrix 
elements of nonlocal operators in the continuum. These nonperturbative matrix elements do not 
suffer from power-divergent mixing on the lattice, which significantly complicates calculations of 
quantities such as the moments of parton distribution functions, provided the smearing scale is kept 
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the Wilson coefficients of the standard operator product expansion and requires the construction of 
a suitable formalism. We demonstrate the feasibility of our approach with examples in real scalar 
field theory. 
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I. INTRODUCTION 

Quantum chromodynamics (QCD) connects hadronic 
matter to its fundamental constituents, quarks, and glu¬ 
ons. Many aspects of QCD are poorly understood, in 
spite of four decades of intense experimental and theoret¬ 
ical effort. As part of this effort, deep inelastic scattering 
(DIS) of leptons from nucleons has played a central role 
in establishing QCD as the reigning theory of the strong 
interaction and continues to serve as a mainstay for at¬ 
tempts to unravel QCD’s mysteries (for complete reviews 
see, for example, Refs. [1] and [2]). 

Theoretically, inclusive DIS cross sections are sepa¬ 
rated into leptonic and hadronic tensors, which capture 
the electroweak and strong dynamics of the scattering 
process, respectively. The hadronic tensor factorizes into 
a convolution of infrared-safe perturbative coefficients 
and parton distribution functions (PDFs), which incor¬ 
porate the low-energy QCD physics and therefore must 
be determined nonperturbatively. PDFs are independent 
of the scattering process-but depend on the target nu¬ 
cleon, while the perturbative coefficients are independent 
of the external scattering states. 

PDFs are important for two reasons. First, they fur¬ 
nish direct knowledge of the constituents of hadronic 
states—the dominant form of visible matter in the 
universe—in terms of the fundamental theory of the 
strong force. Second, they provide constraints on 
hadronic backgrounds at collider experiments, such as 
the Large Hadron Collider. These backgrounds affect 
the sensitivity of a variety of high energy experiments, 
including studies of the properties of the Higgs boson [3- 
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5], and searches for heavy W' and Z' boson production 

[ 6 ], 

PDFs cannot currently be calculated directly from 
QCD using ab initio lattice QCD, because they are de¬ 
fined in terms of light-cone matrix elements that are not 
directly accessible on Euclidean lattices. Hence, PDFs 
are usually extracted from global analyses of different 
experiments [7 15]. 

Naturally a direct nonperturbative computation of 
PDFs is desirable: with sufficient precision, such a calcu¬ 
lation would further constrain global fits in regions that 
may be experimentally inaccessible. 

Until recently, lattice calculations have focused on de¬ 
termining the Mellin moments of PDFs, which are re¬ 
lated to matrix elements of twist-2 operators that can 
be determined on the lattice [16-18]. The lattice regu¬ 
lator breaks Lorentz symmetry, which induces radiative 
mixing between operators of different mass dimensions. 
The resulting power divergences in the lattice spacing 
completely obscure the continuum limit. Moments up to 
the fourth moment can be extracted by carefully choos¬ 
ing the external momenta and operators, but this signif¬ 
icantly limits the precision with which one can extract 
meaningful results for PDFs [19, 20]. Beyond the fourth 
moment, power divergent mixing is inevitable and this 
method cannot provide reliable calculations. 

Ji recently proposed a new approach to directly ex¬ 
tract PDFs from lattice calculations [21, 22], based on 
a large-momentum effective theory [23]. Preliminary re¬ 
sults first appeared in Ref. [24], but a number of issues 
remain, including a complete understanding of the renor¬ 
malization of the relevant lattice matrix elements and the 
practical ability to resolve sufficiently large momenta on 
the lattice. 

Here we propose a new formalism that removes mixing 
in the continuum limit and, in principle, enables the ex¬ 
traction of higher moments of PDFs from lattice QCD. 
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We call this formalism the “smeared operator product ex¬ 
pansion” (sOPE). We expand continuum matrix elements 
in a basis of locally-smeared lattice degrees of freedom. 
The resulting matrix elements are functions of two scales, 
the smearing scale and the renormalization scale. The 
sOPE provides the framework necessary to relate these 
matrix elements, via Wilson coefficients, to phenomeno¬ 
logically useful quantities, such as the hadronic tensor of 
DIS. 

Smearing has been widely used in lattice QCD to re¬ 
duce ultraviolet fluctuations, partially restore rotational 
symmetry and thereby systematically improve the pre¬ 
cision of lattice calculations [25-28]. For a pedagogical 
overview of smearing in lattice calculations, see, for ex¬ 
ample, Ref. [29]. In the sOPE, we implement smearing 
via the gradient flow, a classical evolution of the fields 
in a new dimension, the flow time, toward the station¬ 
ary points of the action [30-34]. Matrix elements deter¬ 
mined nonperturbatively on the lattice require no further 
renormalization, up to a fermionic wave function renor¬ 
malization [30], and remain finite, provided the smearing 
scale is kept fixed in physical units in the continuum limit 
[31, 35]. There are two further advantages: the gradient 
flow allows one to use smearing lengths of only one or 
two lattice spacings, much smaller than hadronic length 
scales on typical lattices, and therefore does not distort 
the low energy physics [36], and also the gradient flow is 
computationally very cheap. 

The gradient flow is now well-established as a tool 
to study lattice gauge theories, with applications from 
scale-setting, i. e. determining the lattice spacing in phys¬ 
ical units [37-41], to studying renormalization in lattice 
gauge theories. For example, the gradient flow has been 
used to define finite-volume renormalization schemes for 
the strong coupling constant [42-46], for operator renor¬ 
malization [47], and to understand the nonperturbative 
scale dependence of renormalized matrix elements [48]. 
Related work has used Ward identities to study chiral 
fermions on the lattice [49] and the energy-momentum 
tensor [50-52], 

In this paper we introduce the gradient flow for a sin¬ 
gle real scalar field with quartic interactions and outline 
some of the properties of the sOPE applied to real scalar 
held theory. Scalar held theories describe some important 
critical phenomena in nature, such as the antiferromag¬ 
netic phase transition, and have a long history as a test¬ 
ing ground for fundamental ideas in quantum held theory 
in four dimensions. For our purposes, the chief advantage 
is the simplicity of Euclidean scalar held theory, which 
lays bare the structure of the sOPE and highlights the 
points of similarity and contrast with the local operator 
product expansion (OPE). 

In the next section, we review Wilson’s OPE and in¬ 
troduce the sOPE. We then apply the sOPE to scalar 
held theory in Sec. Ill, illustrate how the gradient how 
removes power-divergent mixing in Sec. IV, and calcu¬ 
late the perturbative Wilson coefficients to two loops in 
Sec. V. Finally, in Sec. VI, we study the scale dependence 


of the sOPE through renormalization group equations for 
the Wilson coefficients. 


II. THE OPERATOR PRODUCT EXPANSION 

Wilson’s approach to the OPE for a nonlocal operator 
is widely known—see, for example, Ref. [53]—and here 
we review some notation necessary for our discussion. 

We write the OPE for a nonlocal operator, Q(x), as 

Q(x) x ~° 5>Or, M )[0( fe ) ( O)] R . (1) 

k 

The Wilson coefficients c^{x, /i) are complex functions of 
the spacetime separation, x , and renormalization scale, 
/.t, that capture the short-distance physics associated with 
the renormalized local operator [0V)(O)]r. We represent 
renormalized operators by [.. .]r and suppress their de¬ 
pendence on the renormalization scale, /i, for notational 
simplicity. The free-held mass dimension of the local op¬ 
erator determines the leading spacetime dependence of 
the corresponding Wilson coefficient. 

As an example, let us consider the time-ordered two- 
point function of two scalar helds with spacetime sepa¬ 
ration x: T{(j){x)<j)(0)}. In free scalar held theory, the 
OPE is a Laurent expansion. Interactions generate sub¬ 
leading dependence on the spacetime separation in the 
Wilson coefficients, which become functions of the space- 
time separation, the (renormalized) mass Tor, and the 
renormalization scale, /i: 

T{<j>{x)<t>{ 0)} = ^ + c 4> 2 (V x > torx) [cj> 2 (0)]r 

+ 0(x). (2) 

Here we have factored out the leading spacetime depen¬ 
dence from the Wilson coefficients. The 0(x) indicates 
terms of order x, up to logarithmic corrections generated 
by interactions. 

We propose replacing the set of local operators in the 
OPE by their locally smeared counterparts 

Q(x) a ~° ^4(t,;t,/z)S (A:) (t, 0). (3) 

k 

The Wilson coefficients, dk{T,x,n), and the smeared op¬ 
erators, 5^ fc ^(r, 0), are now functions of an extra scale, 
the how time, r. Nevertheless, the leading spacetime 
dependence of the Wilson coefficients is dictated by the 
canonical mass dimension of the corresponding smeared 
operator and is therefore just that of the standard OPE. 

Just as the OPE is only valid for small spacetime sepa¬ 
rations, the sOPE requires small how times. We will see 
that we require that r oc x 2 , where x is assumed to be 
small, to ensure that the Wilson coefficients are indepen¬ 
dent of the external states. 
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For example, returning to the time-ordered two-point 
function, the sOPE is 

T{<t>(x)<t>{ 0)} =- ^- 2 -1 

+ d p 2 (/j,x,[i 2 T, m R x)p 2 (r , 0 ) + 0(x,t), 

( 4 ) 

where we denote smeared fields at flow time r by p(r, a;). 

In the following sections, we will observe four features 
of the smeared Wilson coefficients in the sOPE. 

1. The logarithmic spacetime dependence of the orig¬ 
inal OPE is preserved in the sOPE. 

2. The flow time serves as an ultraviolet regulator for 
the smeared Wilson coefficients, to leading-order in 
perturbation theory. 

3. Beyond leading-order the flow time cannot regular¬ 
ize the Wilson coefficients, because the flow evolu¬ 
tion is classical. We can, however, absorb renormal¬ 
ization scale dependence into the renormalization 
parameters of the original theory. 

4. For any OPE to be meaningful, the Wilson coeffi¬ 
cients must be independent of the external states. 
We can ensure this for the sOPE by choosing s rms < 
X', i.e. the mean smearing radius must be smaller 
than the spacetime extent of the nonlocal operator. 
This ensures that the sOPE remains an expansion 
in approximately local operators. In other words, if 
the gradient flow probes length scales on the order 
of the nonlocal operator, then the sOPE becomes 
a poor expansion for the original operator. 

Although we refer to this expansion as the smeared 
OPE, we really have in mind that the smearing is im¬ 
plemented via the gradient flow. The gradient flow acts 
as a smoothing operation that drives the degrees of free¬ 
dom of the theory to the stationary points of the ac¬ 
tion. In QCD, the gradient flow corresponds to a contin¬ 
uous stout-smearing procedure, an analytic method for 
constructing lattice gauge fields with damped ultraviolet 
fluctuations [25]. The direct comparison of the gradient 
flow to other smearing schemes was first undertaken in 
Ref. [54], 

Many studies of the gradient flow have incorporated a 
small flow-time expansion of fields at nonvanishing flow 
time in terms of local fields at zero flow time [48-52, 55]. 
We can view such an expansion as an OPE in the flow 
time and thereby relate renormalized quantities calcu¬ 
lated at nonzero flow time to the corresponding quanti¬ 
ties in the original theory at vanishing flow time, which 
would otherwise be difficult to compute. 

In this work, we take a different approach. We do 
not expand the flowed fields in terms of original fields at 
nonzero flow time, but rather take as the fundamental 
objects of study the (matrix elements of) fields at pos¬ 
itive flow time. Both approaches reflect the physically- 
motivated expectation that smearing scales much smaller 


than the physical scales of the system should not distort 
the physics in question. The small flow-time expansion 
quantifies any deviations from this expectation and, fur¬ 
thermore, decouples analytic calculations of Wilson co¬ 
efficients in the continuum from lattice calculations with 
smeared degrees of freedom. The sOPE incorporates the 
smearing scale as an inherent scale of the system, which 
requires new Wilson coefficients to be determined. Thus, 
both the small flow-time expansion and the sOPE are 
shaped by the role of the smearing scale as ultraviolet 
regulator and are related, but distinct, conceptual ap¬ 
proaches to the same physics: partially restoring rota¬ 
tional symmetry on the lattice. 

Smearing, in general, is a tool that partially restores 
rotational symmetry [36] and thereby improves the con¬ 
tinuum limit of lattice calculations. Smearing via the 
gradient flow has a number of advantages. In partic¬ 
ular, matrix elements determined nonperturbatively on 
the lattice using smeared degrees of freedom require no 
further renormalization [32], up to fermionic wave func¬ 
tion renormalization [30], and remain finite, provided the 
smearing scale is kept fixed in physical units as the con¬ 
tinuum limit is taken. 

We now turn to a more complete study of the sOPE 
applied to </> 4 scalar field theory, a particularly straight¬ 
forward theory in which to examine the sOPE. We can 
solve the flow-time equations exactly, because the flow¬ 
time evolution is linear in the scalar field. We do not need 
to consider the complications associated with gauge fix¬ 
ing [31, 32]—nor the extra renormalization of fermions 
[30]. Furthermore, for scalar fields the sOPE can be 
understood to be simply a resummation of the original 
OPE. Although it is not necessary for our work, it is also 
interesting to note that the OPE is known to converge for 
Euclidean ^) 4 theory in four dimensions, at a fixed, but 
arbitrary order in the perturbative expansion [56]. This 
result holds at arbitrary spacetime separations, provided 
the external states are of compact support. 

Looking toward future calculations in QCD, we an¬ 
ticipate that the technical issues associated with a flow 
equation that incorporates gauge field interactions will 
slightly complicate the perturbative calculations by in¬ 
creasing the number of diagrams at a given order in 
perturbation theory, but will not alter our conclusions. 
Within the gauge sector, there are no loops of flowed 
fields, because renormalized correlation functions remain 
finite at nonzero flow time [31]. Therefore, at leading- 
order in perturbation theory, the flow time regulates 
ultraviolet divergences; beyond leading-order an appro¬ 
priate renormalization procedure must be incorporated. 
With fermionic fields there is an extra fermion renormal¬ 
ization parameter, calculated in Ref. [30], but this can be 
removed by considering renormalization group invariant 
quantities [30]. We also note that all perturbative calcu¬ 
lations for the sOPE can be carried out in the continuum 
and do not require lattice perturbation theory, which is 
generally more involved [57, 58]. 
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III. GRADIENT FLOW FOR SCALAR FIELD 
THEORY 


On the lattice, we replace the partial derivatives with 
discrete difference operators 


We work in four-dimensional Euclidean scalar field the¬ 
ory with quartic self-interactions and bare mass m, de¬ 
fined by the action 


SM = \f d4 


+ m 2 </> 2 + 


( 5 ) 


To study the sOPE, we introduce the scalar gradient 
flow, which we define through the flow-time evolution 
equation 


dp(j, x) 
dr 


( 6 ) 


where p is a scalar field at nonzero flow time, r, and d 2 
is the Euclidean, four-dimensional Laplacian operator. 
Note that the flow time has units [r] = \x] 2 . 

Imposing the Dirichlet boundary condition p(0, x) = 
4>(x), we may write the exact solution of the flow time 
equation as 

p(r,x) = e Td2 (j)(x), (7) 


or, in the momentum representation, as 

p( T ,p) — e ~ Tp2 4>{p). (8) 


The full solution is 


T lT-un 0 s ) = i • ■ • A /u>0) - traces. (12) 

The spacetime point x is now a node in the lattice, 
x p = an ft, where a is the lattice spacing and is a 
four-component lattice vector. The discrete difference 
operators can be improved to remove discretization ef¬ 
fects, but the simplest such symmetric operator acts on 
a scalar field as 

^ Wx + P) ~ Hx ~ p)), (13) 

where p is the unit vector in the p, th direction. 

The lattice regulator breaks rotational symmetry, 
which induces mixing between twist-2 operators of dif¬ 
ferent mass dimension. On dimensional grounds, the 
mixing coefficients scale with inverse powers of the lat¬ 
tice spacing and these coefficients diverge in the con¬ 
tinuum limit. This is the problem of power-divergent 
mixing on the lattice. For example, the simple bilin¬ 
ear T utt (x) = cj)(x)(l)(x ) mixes with the operator T la ^{x) 
with a coefficient that scales as 1/a 2 . More generally, the 
mixing between T latt ( x ) and any twist-2 operator with an 
even number of derivatives, (a:), scales as l/a 2n . 

The smeared counterparts of these operators, which 
are given by 


P(r,x)= j d 4 y J ^^e ip ' (x v) e Tp2 &(y) 

= 16^2 j d iye ~ (x ~ v)2/4T H y ), ( 9 ) 

which demonstrates explicitly the “smearing” effect of 
the gradient flow: the flow time exponentially suppresses 
ultraviolet modes. We parametrize the smearing radius 
via the root-mean-square smearing length, s rms = V&T. 

The (Euclidean) smeared scalar propagator, for two 
fields at flow times T\ and 72 , is given by 

__ e ~{ ri + T 2 ) k 2 

G p (T U T2,k)= fc2+m2 ■ ( 10 ) 

The flow-time evolution is classical, so any interactions 
must occur at zero flow time. The corresponding Feyn¬ 
man rule for the four-point vertex is just that of the stan¬ 
dard (Euclidean) four-point vertex, V ^ = —A/24. 


(x) = p{r , x)A fil ... A P „p(r, x) - traces. (14) 

do not suffer from this problem. The mixing coefficient 
between £ latt (a;) = p(t,x)p(t,x) and S la4t (x) scales 
as l/r ra , where r is the flow time in physical units. Pro¬ 
vided we keep the dimensionful scale t = ra 2 fixed, where 
t is dimensionless, then as the lattice spacing a decreases, 
the mixing coefficient remains finite, because matrix ele¬ 
ments at nonzero flow time cannot contain any additional 
divergences [31, 35]. 

As a simple illustration of this behavior, let us consider 
the matrix element of two twist-2 operators of different 
mass dimension: 

M cont = <n|0(O)0(O) </>(0)d M cU(0)| Cl). (15) 

This matrix element vanishes for massless scalar fields. 
On the lattice, however, the corresponding matrix ele¬ 
ment, 


IV. MIXING ON THE LATTICE 

Before we examine the sOPE in detail, we demonstrate 
how the gradient flow removes power-divergent mixing on 
the lattice. We consider the example of twist-2 operators 
in scalar field theory, which are symmetric and traceless 
and given by 

T /Jl ...p n (x) = (j){x)dp 1 . ..dfj n <j)(x) - traces. (11) 


Af utt = <ft|0(O)0(O) 0(O)A M A v 0(O)|n), (16) 

is nonzero. 

We show the Feynman diagram for the leading contri¬ 
bution to this matrix element in the left-hand diagram 
of Fig. 1, which is given by 

M (o) = 1 r /a d 4 (afc) kX 

a 4 J-n/a (2^) 4 (P + m 2 ) 2 ’ 


(17) 
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FIG. 1. Diagrams representing the leading contributions to 
the mixing matrix elements with unsmeared fields, , 

and smeared fields, Af/attM- The black solid lines in the 
left-hand diagram represent unsmeared propagators and the 
black solid square and diamond are the operators 0 2 (0) and 
0(0) A m A„ 0(0). The blue dashed lines in the right-hand di¬ 
agram are smeared propagators and the two different blue 
blobs, filled and unfilled, represent the smeared operators 
p 2 (r, 0) and p(r, 0)A^ A„p(r, 0). 


where = (2/a) sin(afc /i /2). We have included a bare 
mass mo to remove spurious infrared divergences, but 
the integral is finite and we can take the massless limit, 
too —t- 0, to match to the continuum massless theory. We 
expand this in powers of the lattice spacing to obtain 


Matt = 


e pir/a j4l, h.2 

_ 1 _ +0(a 2 ) 

4 J_ n/a (2Tr)*(k* + m 2 )* + [ > 

6, 


J [IV 

48M 


k / a 

O(a 0 ) 


(18) 


in the massless limit. This result signals the appear¬ 
ance of power-divergent mixing: the 1/a 2 term diverges 
in the continuum limit. Although this calculation is only 
leading-order in perturbation theory, higher-order contri¬ 
butions do not modify this power-divergent dependence 
on the lattice spacing [36]. 

If we calculate this matrix element with smeared de¬ 
grees of freedom (which we depict in the right-hand dia¬ 
gram of Fig. 1), however, we find 


M, 


(o) 


iatt( T ) = (^Ip( t > 0)p(r, 0) p(r, 0)A Al A^p(r, 0)|f2) 
S uv r /a d 4 fc fc 2 e- 4rfe2 


i 

4 J 

47r 


= 6 


—ir /a (2tt) 4 (P + m 2 )2 

°-T/a 2 _ X 


0{a 2 


fllS 


2567r 2 r 


O(a 0 ). 


(19) 


Once again we have taken the massless limit. In the 
continuum limit, keeping the flow time r fixed in physical 
units, this matrix element tends to a constant, signaling 
the suppression of power-divergent mixing for smeared 
degrees of freedom. 



FIG. 2. Diagrams representing the contributions to the Wil¬ 
son coefficient di at leading-order (left two diagrams) and 
next-to-leading-order (right two diagrams). Solid black and 
dashed blue lines are propagators at vanishing and nonvanish¬ 
ing flow times, respectively. The black squares are unsmeared 
fields 0(0), black dots are interaction vertices at vanishing 
flow time and the blue blob represents the smeared operator 
p 2 (r,0). 


for the leading connected and disconnected operators, 
starting with the disconnected contribution. 


A. Disconnected contributions 


1. Leading order 


We illustrate the leading-order (tree-level) and next- 
to-leading-order (one-loop) contributions to the smeared 
Wilson coefficient for the disconnected operator, I, in 
Fig. 2. As we will see, the leading-order contribution is 
independent of the renormalization scale, /i, because at 
this order the flow time serves as the ultraviolet regula¬ 
tor. Beyond leading order, however, the smeared Wilson 
coefficient will have a renormalization scale dependence. 

We extract the disconnected smeared Wilson coeffi¬ 
cient by considering matrix elements of each of the op¬ 
erators in the sOPE in the vacuum, which removes any 
connected contributions. To 0(x) we have 

<fi|T{0(x)0(o)}|fi) = dl{lix ;^ mx) (n\ I |n) 

+ d p 2(nx, n 2 T,mx)(n\[p 2 (T, 0)] R |n) + 0{x), (20) 

Here we have chosen the normalization of dj so that at 
leading order in the free theory dj is unity. We expand 
each quantity in this expression to one loop according to 


V. WILSON COEFFICIENTS FOR THE SOPE 

The procedure for calculating smeared Wilson coeffi¬ 
cients parallels that for the OPE, discussed in, for ex¬ 
ample, Ref. [53]. With the Feynman rules of Sec. Ill in 
hand, the calculation of smeared Wilson coefficients up to 
next-to-leading-order in perturbation theory is straight¬ 
forward. We determine the smeared Wilson coefficients 


/ = / (0) - A/« + 0(A 2 ), (21) 

where / stands for either a matrix element or Wilson 
coefficient. We have chosen the sign of the one-loop con¬ 
tribution to factor out the sign of the coupling constant 
arising from the Feynman rule for the four-point vertex 
F( 4 ) = —A/24. 

The tree-level Wilson coefficient is then given by the 
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small spacetime behavior of 


2. Next-to-leading-order 


=° {<Q|T«*).«0)}|S!> 

-<Sl|^(r,0)|!l>}™ mii , (22) 

where we have neglected the arguments of the coefficient 
for clarity and used the fact that the tree-level connected 
coefficient is djv* = 1. The subscript indicates that we 
must expand the result to 0(m 2 ). For more details about 
the calculation of Wilson coefficients see, for example, 
Ref. [53]. The corresponding Feynman integral represen¬ 
tation is 


di 0) x =° 47 tV 


= 47r 2 a; 2 


^ik-x g—2 k z 


k 2 


Ak-x —2kr 

e — e 


0(m 2 ) 

1 m 2 


, (23) 


where 


d 4 fc 

(2^ 


(24) 


The smeared Wilson coefficient is 


c4 0) = l- 


8 T 


x 2 \ 

7E-l+log( — J 


(25) 


with 7 e — 0.577216 the Euler-Mascheroni constant. 

We can compare this result with the Wilson coefficient 
for the OPE for d = 4 — 2e dimensions [53]: 


.(0) _ i , m2x2 

c i - 1 + , 


- + 1 + 7 E + log 

e 


2 2 
7 TfX X 


• (26) 


At one loop, the smeared Wilson coefficient is given by 
d (1) r 

= {(S!|T{#t)«0)}|S2) 

-<n|p"(r,0)|n)}[V (28) 

The four-point interaction in this diagram, which we 
show in Fig. 2, appears at zero flow time. Therefore the 
flow time cannot act as a regulator for the momentum 
integral over and we must introduce a renormaliza¬ 
tion procedure. We use dimensional regularization and 
the MS scheme. The double integral is straightforward, 
however, because the two integrals can be carried out 
separately: 

4 1} _ r r e iki X - e -2fc i r i r 1 l 

4? t 2 x 2 \J kl (k 2 + m 2 ) 2 2J k2 k% + m 2 J 0(m2) ‘ 

(29) 

We find, for m > 0, 


7(P 


2 2 
m z x z 

1287T 2 


7 E - 1 + log 



1 + log 



Here we see that the second term, which is a function 
of the new renormalization scale, /i, and the bare mass, 
?n, is nothing other than the one-loop contribution to the 
mass renormalization of the original theory, in the MS 
scheme [59]. Moreover, the factor containing the flow 
time is identical to the 0{m 2 ) term from the leading- 
order contribution, d[°\ in Eq. (25). 

Therefore, we can simply combine the leading-order 
and next-to-leading-order terms to give 


2 2 2 
, „ ar mtx z 

dl = 1 -8 7 + ^ 


7 E - 1 + log 



0(A 2 ), 

(31) 


Then in the MS scheme, this becomes 


( 0 ) _ 1 rn 2 x 2 

c i - 1 + 4 


1 + 27 E + log 


2 2 
fL X Z 

16 


(27) 


Although the finite contribution to these expressions 
cannot be directly compared, because we have expressed 
the Wilson coefficients in two different renormalization 
schemes, we note three important features. First, the 
logarithmic dependence on the spacetime separation is 
identical. Second, the flow time r plays the role of the 
renormalization scale at leading order. Third, we see 
that for small spacetime separations, we require a small 
flow-time parameter. If we do not choose the flow-time 
parameter appropriately, we generate large logarithmic 
contributions to the smeared Wilson coefficients and the 
sOPE exhibits poor convergence properties, even at small 
spacetime separations. 


where tor is the renormalized mass in the MS scheme, 
given by = .21“ ^n 2 and Z m is the mass renormaliza¬ 
tion [59]: 


Z m — 1 + 


A 

167T 2 


1 + log 



+ 0( A 2 ). 


(32) 


This is a clear, next-to-leading-order example of how 
the divergences of the theory at nonzero flow time are ab¬ 
sorbed by the renormalization parameters of the original 
theory at zero flow time. In other words, the renormal¬ 
ized theory at zero flow time remains ultraviolet finite at 
nonzero flow time. 


B. Connected contribution 

We illustrate the leading- and next-to-leading-order 
contributions to the connected Wilson coefficient 



7 



FIG. 3. Diagrams representing the leading-order (one-loop) 
contributions to the Wilson coefficient d p 2 . Details of the 
Feynman diagrams provided in the caption of Fig. 2. 



FIG. 4. Diagrams representing the contributions to the next- 
to-leading-order (two-loop) contributions to the Wilson coef¬ 
ficients d p 2 . For details of the Feynman diagrams, see the 
caption of Fig. 2. 


d p 2 ( fix , /x 2 r, mx) of Eq. (4) in Figs. 3 and 4, respectively. 

Throughout this section, we neglect diagrams that are 
trivially incorporated as part of the wave function renor¬ 
malization of the external fields, such as the one-loop ex¬ 
amples illustrated in Fig. 5. Provided the original theory 
at zero flow time is renormalized, counterterms cancel 
these contributions completely. 


1. Leading order 


We extract the leading-order connected Wilson coeffi¬ 
cient by considering matrix elements of each of the oper¬ 
ators in the sOPE coupled to two external fields, which 
removes any disconnected contributions. We can then 
read off the one-loop contribution to the Wilson coeffi¬ 
cient, </ 2 ^, by matching terms at 0(A). This contribution 



FIG. 5. Example diagrams that are naturally incorporated 
in the renormalized external propagators. We do not include 
these contributions in our explicit calculations of the Wilson 
coefficients. For details of the Feynman diagrams, see the 
caption of Fig. 2. 


is given by the small spacetime behavior of 

x =° {(n\T{<i>(x)m}t(pS(P2m) (1) 

- (n|p 2 (r,0)^(pi)^(p 2 )|fi) (1) } ■ (33) 

The corresponding Feynman integral is 

1 f 1 f e ik ' x — e - ( fe2+<?2 ) r 1 

(p 2 + m 2 )(p 2 + m 2 ) | 2 J k (k 2 + m 2 )(q 2 + m 2 ) J ’ 

(34) 


where q = k — p± — p 2 . We extract the smeared Wilson 
coefficient by examining the small spacetime behavior of 
the integral in curly braces, expanded to O(m 0 ): 


d 


(1) x~0 
nl ~ 


1 r gik-x _ e ~(k 2 +q 2 )r ) 

2 J k (.k 2 + m 2 )(q 2 + to 2 ) J 

J (9(m u 


(35) 


The smeared Wilson coefficients must be independent 
of the external states to ensure that the sOPE is truly an 
operator expansion. In this particular case, we require 
that d^J is independent of the external momenta p\ and 
P 2 and the mass. By taking a derivative with respect to 
one of the external momenta, 


d (i) _ f Qi [e lk ' x - e ( fe2 +9 2 ) T (l + (g 2 + to 2 )t)] 
dp* p2 Jk (fc 2 + m 2 )(q 2 + to 2 ) 2 

(36) 

we obtain a convergent integral. The x —> 0 limit of this 
integral is now well defined and only vanishes if the flow 
time is related to the spacetime separation. An analo¬ 
gous result holds if we take a derivative with respect to 
the external mass, to. On dimensional grounds, then, we 
guarantee that the smeared Wilson coefficient is indepen¬ 
dent of the external states by choosing r oc x 2 . 

On physical grounds, however, we require that the 
smearing radius is smaller than the spacetime extent of 
the nonlocal operator: s rms < x. This choice ensures that 
the sOPE remains an expansion in local operators. Phys¬ 
ically speaking, if the gradient flow probes length scales 
on the order of the nonlocal operator, which we represent 
in Fig. 6, then smeared operators would cease to be (ap¬ 
proximately) local. In other words, the sOPE would be 
a poor expansion for the original operator. This is the 
physical origin of the third feature that we saw in the 
leading-order disconnected contribution, Eq. (25): small 
spacetime separations require small flow times to ensure 
convergence. The OPE requires small spacetime separa¬ 
tions for good convergence, so it follows that the sOPE 
requires small flow times as well. 

There is a complementary viewpoint that elucidates 
the small flow-time requirement more quantitatively: the 
small flow-time expansion [31, 50, 52]. In this case we 
see that the derivative of Eq. (36) is independent of the 
external momenta in the x -4 0 limit, up to terms linear 
in the flow time: 


_d_ 

d pi 



qi(e ik ' x — l) 


Ik (^ 2 + rn 2 ){q 2 + wi 2 )" 


+ 0(t). 


(37) 










FIG. 6. At small flow times (left-hand diagram), the smeared 
operators are localized relative to the spacetime separation of 
the nonlocal operator. At large flow-time values (right-hand 
diagram), however, the smearing radius probes the scale of the 
nonlocal and the sOPE is a poor approximation to the original 
nonlocal operator. For details of the Feynman diagrams, see 
the caption of Fig. 2. 


From this, it is clear that the Wilson coefficients will 
be independent of the external states, up to terms lin¬ 
ear in the flow time. Therefore, provided the flow time 
is small relative to the spacetime separation, the sOPE 
is an expansion in approximately local operators, in a 
quantifiable sense. Moreover, the flow-time dependence 
will cancel in the product of the Wilson coefficient and 
its associated matrix element, to the desired order in the 
flow time. 

We are free to set p± = P 2 = 0 in Eq. (35), because the 
smeared Wilson coefficient is independent of the external 
momenta and the mass to the order at which we are 
working [53]. Expanding in the mass, we obtain 



FIG. 7. Two-loop diagrams that appear at 0(m 2 ) and there¬ 
fore do not contribute to d p 2 . For details of the Feynman 
diagrams, see the caption of Fig. 2. 


is 


(fi|p 2 (r,0)</>(pi)(Xp2)|ft) = 1 

+ 3^ t 1 +7E + log ( 2m2r )] + 0 (t,X 2 ) (41) 

for sufficiently small flow times. Here we have dropped 
the external fields, which we take to be on shell, for sim¬ 
plicity. The product of this matrix element with the Wil¬ 
son coefficient is independent of the flow time to one loop 
and 0(t), as we would expect: 


d p 2(H|p 2 (r,0)^(pi)</)(p2)|^) = 1 

A 


+ 


32t r 2 


2q E + log 


+ 0 (t, A 2 ). (42) 


2. Next-to-leading-order 


d« = i I 6 '"- 6 


-2 k A 


32t r 2 


(fc 2 ) 2 

TE-i+togfe) 


(38) 


Combining this with the leading-order contribution, 
which is just unity, we have 


dp2 ~ 1+ 32^ 


7 E - 1 + log 


_2 \ "I 


8 T 


+ 0 ( A 2 ). (39) 


In contrast, the Wilson coefficient for the OPE in the 
MS scheme, denoted by c, is 


To determine the next-to-leading-order contribution to 
the connected Wilson coefficient, we must incorporate the 
Feynman diagrams of Fig. 4. The two-loop diagrams of 
Fig. 7 do not contribute to the Wilson coefficient d p 2 , 
because they appear at (2 (to 2 ). 

Looking at the Feynman diagrams of Fig. 4, and bear¬ 
ing in mind our experience of the next-to-leading-order 
contribution to di, we immediately observe that these di¬ 
agrams are simply the product of the one-loop Wilson co¬ 
efficient, d^V, and the next-to-leading-order renormalized 
four-point vertex. This contribution to the vertex is, of 
course, nothing other than the next-to-leading-order con¬ 
tribution to the renormalized coupling constant. Thus, 
we write the Wilson coefficient quite simply as 


f-;iy 


= 1 


32?r 2 


1 + 27 E + log 


2 2 
[TX* 

16 


0( A 2 ). (40) 


dpi — 1 + 


Ar 

32tt 2 


7 e - 1 + log 



+ 0(A 3 ), 


(43) 


We note the occurrence of the three features we observed 
for the leading-order disconnected contribution: the same 
spacetime dependence for both smeared and unsmeared 
coefficients; the appearance of the flow time as a leading- 
order regulator; and the need to choose a small flow time 
for small spacetime separations to avoid large logarith¬ 
mic contributions. From the small flow-time expansion 
viewpoint, we can confirm that the flow-time dependence 
ultimately cancels to the desired order. For example, the 
matrix element of p 2 (r, 0) coupled to two external fields 


where Ar is the renormalized coupling. In the MS 
scheme, the renormalized coupling is given by 

A » = A -2(I^ tog (£) +0 ( AS >’ <«> 

For a calculation of this to five loops, see [59]. 

As we move beyond the leading-order Wilson coeffi¬ 
cients, i.e. dj 0 ' 1 and c/V, divergent radiative corrections 
appear in our calculations, because all field interaction 
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vertices appear at zero flow time. These divergences 
can be removed by the renormalization parameters of the 
original theory and the renormalization scale dependence 
of the smeared operators is completely contained in the 
renormalization parameters of the original theory. This 
is to be expected: it follows from the fact that renormal¬ 
ized matrix elements at zero flow time remain finite at 
nonzero flow time and require no further renormalization 
[31, 35], 

Ultimately, for realistic calculations in lattice QCD, 
the perturbative calculation of smeared Wilson coeffi¬ 
cients must be combined with nonperturbative determi¬ 
nations of matrix elements at hadronic energy scales. 
Scalar </> 4 theory is not asymptotically free in four dimen¬ 
sions, but we can understand the mathematical features 
of the sOPE in more detail by studying the renormal¬ 
ization group equations for the simple example of scalar 
fields. 


VI. RENORMALIZATION GROUP EQUATIONS 


We consider the matrix elements of scalar operators 
coupled to N external, unsmeared scalar fields. We can 
derive renormalization group equations for the Wilson 
coefficients of the OPE by considering the scale depen¬ 
dence of suitably chosen Green functions [53]. The Green 
function for N + 2 external scalar fields obeys the renor¬ 
malization group equation 

+ (iV + 2)7 

dp 

while the Green function of the renormalized operator 
(/> R (0) coupled to N external scalar fields satisfies 


(fi|<A(Pi)---0(pjv+2)|fi> =0, (45) 


p-— 7 m + iVq 

d fi 


(Q|^(o)0(pi)... <j){p N )\n) = o. 


(46) 

Here the renormalization group operator for scalar field 
theory is 


d d 


d 

% = ^ 

\,m R 

7mWR a 

fi,m R 


(47) 




and the coefficients are 


/3 ms (A) 

7^(A) 

7 M (A) 


dA _ 3A 2 

^dp (167T 2 ) 2 
p d log(m R ) 


0 ( A 3 ), 


2 dp 
A 

^ 2(16tt 2 ) 
pdlog (Z#) 
2 dp 


5A 2 

12(16tt 2 ) 2 

A 2 

12(16?r 2 ) 2 


fO(A 3 ), 
+ 0(A 3 ). 


(48) 

(49) 

(50) 


In general these coefficients depend on the mass, the 
renormalization scale, and the renormalized coupling 


constant, but in a mass independent renormalization 
scheme, such as the MS scheme, they depend on the 
renormalization scale only through the renormalized cou¬ 
pling constant. These functions are known to five loops 
for the (J(IV)-symmetric theory, given by the V-multiplet 
$ = {(/>!,...,<M [59, 60], 

Returning again to our example of the OPE for the 
two-point function, Eq. (2), we can derive an renor¬ 
malization group equation for the Wilson coefficient 2 
by coupling these operators to two external scalar Helds 
and using Eqs. (45) and (46) (for further details, see, for 
example, Ref. [53]): 

d , , 

^d/i +2 W + 7m) 

Just as we might expect, the anomalous dimension of 
the Wilson coefficient 2 is equal to the difference be¬ 
tween the anomalous dimension of (f>(x)(j)(0) and that of 
[^ 2 (0)]r. 



A. Renormalization group equations for the sOPE 
Wilson coefficients 


The flow time, r, introduces a new scale into the prob¬ 
lem. In principle one can view the flow time as just 
another external scale, unrelated to the renormalization 
scale p. In this case, it is natural to modify the renor¬ 
malization group equations to account for the change in 
the Green functions as we change both p and r: 


d d d 

p-- > p--2r —. 

d fi d fi dr 


(52) 


Our choice of differential operator is constrained by 
the mass dimension of each scale, p and r, but is not 
unique. In particular, one could choose the operator 
pd/dp + rd/dr. This freedom does not affect the logic 
of our discussion nor our conclusion, because alternative 
conventions can be absorbed into the definition of the 
renormalization parameters and anomalous dimensions 
in, for example, Eq. (58). 

At this stage it is worth commenting on the two scales 
in the problem, p and r. In DIS, the spacetime separation 
of the corresponding OPE, x , and renormalization scale 
p are, in principle, two distinct scales. The spacetime 
separation is provided by the inverse momentum transfer 
of a particular DIS experiment or set of experiments. The 
renormalization scale, however, is a theoretical choice, 
and ultimately physical quantities should not depend on 
the renormalization scale. It is generally convenient to 
choose p = 1/a;, but it is not strictly necessary. 

The relationship between the flow time and the renor¬ 
malization scale is analogous and these two scales are 
distinct. For the sOPE, the flow time can be considered 
as simply an external scale, imposed by some particular 
lattice “experiment”, and the renormalization scale is a 
convenient theoretical choice. We will see that it is help¬ 
ful to tie these scales together, to reduce the two-scale 
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problem to a single scale, but this is not formally neces¬ 
sary. Therefore, in the following analysis, the flow and 
renormalization scale should be understood as completely 
independent scales. 

Considering again the sOPE for the two-point func¬ 
tion, Eq. (4), we determine the renormalization group 
equation for the smeared Wilson coefficient, d p 2 , by fol¬ 
lowing a procedure analogous to that outlined above. 

We assume that we are working in the small flow-time 
regime, which allows us to relate operators at vanishing 
and nonvanishing flow time [31, 50, 52]: 

[</» 2 (0)] r = Z p 2 0) + 0(t). (53) 

This coefficient satisfies 

log(Z p2 (rV)) = 2y m , (54) 

and to one loop is given by 

Z p 2 (t,h 2 ) = 1 - [1 +7e +log(2r/r 2 )] +C>(A 2 ,t). 

(55) 

We apply the renormalization group operator 

My- - 2t 4~ + (N + 2)y (56) 

ay dr 

to matrix elements of the operators in Eq. (53) coupled 
to N external scalar fields to obtain 

/y- 2T s +2(< >- 7) 

Here 

Cp 2 = T ^ r log {Z P 2 (r, y 2 )) (58) 

is an anomalous dimension associated with the flow-time 
dependence of the operator p 2 (r, 0). The renormalization 
group equation for the corresponding matrix element of 
p 2 (r, 0) coupled to N external fields is given by 

^ 2r^ + 2 Cy + IV 7 

x (^Im 2 (t, 0 )^(pi) . ..(j){p N ) |fi) = 0(t). (59) 

We note that this equation only holds provided the flow 
time is small compared with the momenta of the external 
particles, which for DIS would be of the order of hadronic 
scales. 

If we now demand that the smearing scale, r, and the 
inverse of the renormalization scale, y, are proportional 
to each other, i.e., r = b/y 2 with b real, then the renor¬ 
malization group equation becomes 

2^— + 2£ p 2 + IV 7 

x (^ 1 P 2 (b/y 2 , 0) <Xpi). ..<i>(p N )\n) = 0(b). (60) 


d p 2 = O(r). (57) 


This renormalization group equation provides the start¬ 
ing point for a nonperturbative step-scaling method 
[47, 61] that evolves nonperturbative matrix elements to 
a high scale, where they can be combined with pertur¬ 
bative smeared Wilson coefficients. Here the renormal¬ 
ization group equation holds in the small flow time limit, 
or in other words, provided b <C 1. This constraint auto¬ 
matically ensures that the flow time is also smaller than 
any hadronic length scales, Aqcd <C y 2 <C 1/r, and is 
generally true for practical step-scaling methods [47, 61]. 
We are currently investigating this approach in QCD. 


VII. CONCLUSION 

We have proposed a new method, the smeared oper¬ 
ator product expansion (sOPE), to extract matrix ele¬ 
ments from numerical nonperturbative calculations with¬ 
out power divergent mixing. The smeared operator 
product expansion is a general framework relevant to 
any asymptotically free theory with nonperturbative ma¬ 
trix elements that suffer from power divergent mixing. 
Within QCD, the most obvious application is to deep 
inelastic scattering, but other applications include non¬ 
perturbative determinations of K —> tttt decays [62, 63] 
and H-meson mixing [64]. Beyond QCD, applications in¬ 
clude nonperturbative studies of critical phenomena in 
the Heisenberg model, spin systems and other condensed 
matter systems. 

In the sOPE, we expand nonlocal operators in a basis 
of smeared operators, the matrix elements of which can 
be determined on the lattice. We implement the smear¬ 
ing via the gradient flow, a classical evolution of the the¬ 
ory in a new dimension that smooths ultraviolet fluctu¬ 
ations. The continuum limit of these matrix elements is 
free of power divergent mixing, provided the localization 
scale, the smearing length, is kept fixed in the continuum 
limit. The resulting matrix elements are functions of two 
scales, the renormalization scale and the smearing length. 
The sOPE systematically relates these matrix elements 
to smeared Wilson coefficients, which can be calculated 
in perturbation theory, thereby providing a complete de¬ 
termination of the nonlocal operators. 
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